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\,  SUMMARY 

n( 

The  solution  of  the  subsonic  flutter  problem,  when  the  commonly  used 
linear  differential  equation  model  is  replaced  by  the  more  correct  linear  integro- 
differential  equation  model,  is  studied  and  the  nature  of  the  system's  free 
motion  established.  The  different  forms  appropriate  to  two-dimensional  and  three- 
dimensional  flow,  and  to  the  cases  when  the  system  has  a  zero  characteristic 
value,  are  developed  in  detail.  It  is  shown  that,  for  large  time  t  ,  the 
behaviour  can  variously  be  like  Hriog^t ,'  or  t  3 
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THE  FORM  OF  THE  SOLUTIONS  OF  THE  LINEAR  INTEGRO-DIFFERENTIAL 
EQUATIONS  OF  SUBSONIC  AEROELASTICITY 


0.  L,  Woodcock 


CORRIGENDA 


1  The  analysis  of  this  paper  has  been  based  on  the  assumption  that  the 

Laplace  transform  g(^)  of  that  part,  G(vt),  of  the  indicial  aerodynamic  matrix 

(equation  (2))  which  remains  when  the  constant  and  impulsive  terms  are  excluded, 

2  10 
had  the  form  given  by  Milne  based  on  the  work  of  Garner  and  Milne  .  This  form, 

reproduced  in  equation  (11)  is 


5d)  ■  t 


where  N^  is  zero  except  in  the  two-dimensional  case.  However  a  closer  examina¬ 
tion  of  Refs  2  and  10  reveals  that  their  authors  only  showed  the  existence  in  the 
three-dimensional  case  of  the  constant  and  (£)  log  (£•)  terms  in' (11).  Moreover, 

17  v  v  . 

using  the  series  expansion  for  the  three-dimensional  kernel  function*,  one 
finds  that  G  has  the  form 


*  The  kernel  is  a  power  series  in  p  ,  plus  p  log  p  times  another  power 
series  in  p  . 


(C-l) 


5®  -  E  Z  Z  i°8W,"Z, 

s=0  r=l  s=r 

Further  confirmation  of  the  existence  of  the  higher  powers  of  log(^)  is  obtained 
when  we  expand  the  known  form  (equations  (12)  and  (13))  of  G  (£■)  for  the  two- 
dimensional  incompressible  case.  The  series  then  has  the  form 


m  ■  E  ?)\  *  E  E  •  «-2> 

s=0  r=l  s=0 

1 8 

The  same  form  of  series  is  also  obtained  when  one  uses  the  series  expansion  of 
the  kernel  function  for  two-dimensional  compressible  flow. 

Taking  (C-2)  as  the  general  form  instead  of  (11),  and  noting  that  N 

IT  S 

is  zero  in  the  three-dimensional  case  when  r  >  s  ,  the  asymptotic  expansion  of 
G(vt)  is  changed  from  (3)  to 


G(vx) 


~  ZZg 


logr(vx) 

.  . r+s+1 

(vt) 


(C-3) 


where,  in  the  three-dimensional  case  G  is  zero  for  r  *  s  .  The  coefficients 

rs 

in  (C-3)  are  not  as  simply  related  to  those  of  (C-2)  as  were  those  of  (3)  to 

those  of  (11).  Any  particular  G  is  a  linear  combination  of  the 

rs 

Nr+U+]  g_u  (u  =  0  +  s).  For  example  we  have 


00 

G0> 

G.0 

G.l 

etc. 


-  N 


10 


NH  +  (2  -  2Y)N20 


-  2N 


20 


4N2]  +  (18  -  12y)N3Q 


J 


(C-4) 


This  change  -  (C-2)  instead  of  (11)  -  makes  no  change  to  the  character  of 
the  subsequent  analysis;  and,  in  particular,  no  change  in  every  case  to  the  form 
of  the  dominant  term,  for  large  time,  of  the  fundamental  solution  Xn  .  There  is 


some  apparent  modification  of  the  overall  form  of  Xq  in  the  three-dimensional 
case.  We  find  instead  of  (46) 


v  £  W.  1  +  v3  ^logr(vx) 

i=l  s=3  r-0 


5r,(T> 

(vt)2"8 


(C-5) 


The  change  of  form  is  only  apparent,  for  the  fact  was  overlooked  that  the  S 
of  equation  (46)  were  zero  for  s  <  (r  +  3)  ;  and  so  putting  s  =  (u  +  r)  in 
(46)*  transforms  it  into  a  form  like  (C-5). 


When  the  characteristic  equation  has  a  zero  root  the  form  is  still  that 
given  in  equation  (47)  though  the  coefficients  Urg  ,  similar  to  the  Srg 
appearing  in  (C-5),  are  now  derived  differently.  In  the  two-dimensional  case  the 
fundamental  solutions  are  those  given  in  equations  (39)  and  (51)  with  a  similar 
proviso  about  the  coefficients. 


2  Additional  references 


No.  Author 


Title,  etc 


17 


18 


C. E.  Watkins 
H.L.  Runyan 

D. S.  Wools ton 


On  the  kernel  function  of  the  integral  equation 
relating  the  lift  and  downwash  distributions  of 
oscillating  finite  wings  in  subsonic  flow. 

NACA  Report  1234  (1955) 


Deborah  J.  Salmond  Evaluation  of  two-dimensional  subsonic  oscillatory 

airforce  coefficients  and  loading  distributions. 
RAE  Technical  Report  79096  (1979) 


3  Corrections  to  the  paper  as  written,  %e  without  making  the  change 

suggested  in  section  1  above. 


p5  equation  (10)  -  replace  v  -*■  <*>  by  v  -*■  “ 

2  2 

p7  equation  (16)  -  replace  v  by  v 

p8  line  before  equation  (22)  -  insert  at  end:  ,  with  the  wing  chord  as  the 
reference  length 

line  after  equation  (22)  -  add  asterisk  to  choice  and  insert  footnote: 

*  For  finite  aspect  ratio  wings  appropriate  values  of  pg  will  be 
larger.  The  results  of  Ref  9  suggest  for  example  a  value  roughly 
twice  as  big  for  an  aspect  ratio  4  wing. 

p 1 2  first  line  -  replace  'T_  (s  =£  0)  ’  by  'T  for  r  <  s' 

U8  IT  8 


The  upper  limit  of  the  third  summation  in  equation  (46)  should  be  (s  -  3). 


•  m 


s 


pi 2  equation  (38)  -  replace  m  =  1  by  m  -  1  as  the  upper  limit  of  the  second 
summation 

equation  (39)  -  replace  s  by  00  as  the  upper  limit  of  the  third  summation 

line  after  equation  (42)  -  replace  'S  0  are  zero'  by  'S  are  zero  for 

r  ^  rs 

s  <  (r  +  3)'.  Replace  Tft  by  'T  for  r  <  s' 

us  rs 

equation  (43)  -  replace  r  +  s  -  1  by  s  -  1  as  the  upper  limit  of  the 
summation 

pl3  third  line  after  equation  (45)  -  replace  infinite  by  infinity 

pl4  third  line  of  section  4.2  -  replace  'each  matrix  »'  by  ',  for  s<  (r  +  3), 

each  matrix  S  ' 
rs 

equation  (46)  -  replace  s  by  s  -  3  as  the  upper  limit  of  the  third 
summation 

equation  (47)  -  replace  s  by  “  as  the  upper  limit  of  the  third  summation; 
2 

and  replace  log  (vt)  by  log  (vt) 

pi 5  equation  (51)  -  replace  s  by  ”  as  the  upper  limit  of  the  last  summation 
pl7  line  8  -  replace  (14)  by  (16) 

second  line  after  equation  (59)  -  replace  'in'  by  'is',  and  'had'  by  'hand' 
p29  first  line  -  replace  'Wang'  by  'Wong* 
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1  INTRODUCTION 

Traditionally  the  flutter  problem  has  been  treated  as  the  determination  of 

1-4 

the  stability  of  a  system  of  linear  differential  equations.  Certain  authors 

have  however  been  aware  of  the  fact  that  a  truer  representation  of  the  physical 

system,  for  small  perturbations  from  an  equilibrium  state,  is  a  system  of  linear 

integro-differential  equations,  and  have  made  some  attempt  to  assess  its  signifi- 

13  4 

cance.  One  way  of  doing  this  *  ’  is  to  approximate  to  the  system  of  integro- 
differential  equations  by  a  higher  order  system  of  differential  equations  but  this 
involves  an  approximation  to  the  indicial  aerodynamic  matrix  which  loses  some  of 
its  significant  characteristics. 

Stability,  as  the  flutter  investigator,  or  indeed  the  worker  in  the  field  of 
active  control  or  aircraft  stability,  understands  it,  is  a  property  of  systems  of 
differential  equations.  When  one  comes  to  integro-differential  equations  one  has 
to  redefine  the  concept  of  stability  if  it  is  going  to  be  used.  It  is  possible, 
for  example,  as  we  shall  see  subsequently,  to  have  such  a  system,  whose  response 
to  an  instantaneous  excitation  of  any  sort  always  ultimately  dies  away,  and  which 
yet  will  respond  to  a  certain  non-ins tantaneous  excitation,  which  ultimately 
becomes  zero,  with  a  response  which  does  not  tend  to  zero.  This  is  a  consequence 
of  the  hereditary  nature  of  such  systems.  Do  we  therefore  say  that  such  a  system 
is  asymptotically  stable  if  the  response  to  any  instantaneous  excitation  ultim¬ 
ately  tends  to  zero,  or  do  we  require  there  to  be  no  excitation  ultimately  tending 
to  zero  which  will  give  a  response  which  does  not  ultimately  tend  to  zero?  The 
precise  definition  of  stability  thus  indicates  a  basic  difference  between  the 
differential  equation  and  integro-differential  equation  models.  The  approximation 
of  Refs  1,  3  and  4,  while  providing  a  much  improved  approximation  to  the  indicial 
aerodynamic  matrix  over  that  implicit  in  the  traditional  approach,  does  eliminate 
this  basic  difference. 

The  only  author  to  consider  the  exact  solution  of  the  linear  integro- 

2  6 

differential  equations  of  flutter  theory  is  Milne  *  .  The  present  study  is  a 
contribution  to  the  same  subject  from  a  somewhat  different  viewpoint. 

3 

This  paper  is  written  as  a  companion  to  papers  by  Lawrence  and  Jackson  , 

4  #  9  B 

and  Woodcock  and  Lawrence  .  Originally  it  was  intended  that  these  two  papers 
should  be  combined  and  issued  as  an  ARC  R  &  M.  The  present  note  is  largely  a 
basic  study  which  it  was  planned  also  to  include  in  the  R  &  M. 

2  EQUATIONS  OF  MOTION 

For  small  perturbations  the  free  motion  of  a  deforming  lifting  surface,  in 
subsonic  flow  is  given  by  the  solution  of  the  matrix  integro-differential 
equation^. 


I 


4 


(A  -  Aj)  +  Eq 
dx 


=  V2  j  K 


dq(T0) 

K(T  '  V  “d TT"  dTo 


(0 


where 


(A  -  Aj) 
E 

A1 

KCx) 


q  (t) 


is  the  structural  inertia  matrix 
is  the  structural  stiffness  matrix 
is  the  aerodynamic  inertia  matrix 
is  the  indicial  aerodynamic  matrix 

is  the  column  matrix  of  the  generalised  coordinates. 


Structural  damping  has  been  neglected  but,  if  desired,  one  can  approximately* 
take  account  of  it  by  including  a  term  D  ^  on  the  left  hand  side,  where  D  is 
the  structural  damping  matrix.  The  indicial  aerodynamic  matrix  K(t)  ,  which 
results  from  a  displacement  with  time  dependence  H(t)  ,  has  the  form 


K(t) 


6  (t 

l0 


V*(T) 

2 

V 


(2) 


where  H(t)  ,  6(x)  are  the  right-handed  Heaviside  step  and  Dirac  delta  functions 
(of  Appendix  B) ;  Aq ,  Aj ,  are  constant  matrices;  and  G(vt)  has  an  asymptotic 
expansion  of  the  form**  (af  Refs  2  and  10) 


G(vt) 


00 


(3) 


Thus  is  the  ultimate  value,  for  large  time,  of  the  indicial  aerodynamic 

matrix.  If  q  =  qe1WT  ,  where  w  is  real,  and  t  is  infinitely  large  (ie  simple 
harmonic  motion  of  infinite  duration)  the  coefficient  of  q  on  the  right  hand 
side  of  (1)  becomes  the  oscillatory  aerodynamic  matrix  v^Aj  ^2"~~  b(^)  , 

where  B  and  C  are  real  for  real  argument.  The  Laplace  transform  K(p)  of 
K(t)  is  therefore  related  to  these  oscillatory  aerodynamic  matrices  by 


*  The  term  structural  damping  embraces  a  multitude  of  effects:  material  damping, 
friction,  backlash  etc;  and  a  truer  representation,  if  possible,  would 
introduce  considerable  complication. 

**  Nq  is  zero  except  in  the  two-dimensional  case. 
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pK(p)  =  "  ■|(^}  Aj  +  £  B(-  +  C(-  (4) 

where  B  and  C  are  now  continued  analytically  into  the  complex  plane. 
Relationships  therefore  exist  between  the  matrices  appearing  in  (2)  and  the 
oscillatory  aerodynamic  matrices.  Using  the  theorem  that 

•^tlpyCp)!  =  y(“)  (5) 

p+0*  ’ 

where  y  is  any  function*  and  y  its  Laplace  transform,  we  immediately  see 
that 

Ka  =  -  C(0)  =  (say)  -  CQ  .  (6) 

Now,  from  (2) 

pit(p)  =  k,  .  £{a0  ♦  g(£)}  -  (£)2a,  (7) 

and  so  (4)  can  be  rewritten 


If  now  we  go  to  the  limit  p  =  ioo  we  get  the  relationships 

Aq  +  «{G(i°c)}  =  -  B(°°)  =  (say)  -  B^  (9) 

[vj?{G(iv)}l  =  CQ  -  C(»)  =  (say)  CQ  -  .  (10)' 

From  the  work  of  Milne^  we  know  that  g(£)  taken  as  single  valued  in  the 
complex  plane  cut  along  the  negative  real  axis,  has  the  form  (of  equation  (3)  and 
Appendix  A) 

00  00 

s®  -  I  I  @'\  <"> 

s=0  s=0 


*  Strictly  any  function  which  is  bounded  for  real  t  >  0.  We  apply  the  theorem 
to  the  last  term  of  (2). 
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where  Nq  is  zero  except  in  the  two-dimensional  case.  In  the  particular  case  of 
two-dimensional  incompressible  flow  it  is  easily  seen,  from  the  known  analytical 
solution  ( eg  Ref  12)  that  g(£)  has  the  form 


02) 


where  R,  S,  T  are  constant  matrices  and  2  is  the  Theodorsen  function 


a©  = 


Ko^y +  ki  yy 


(13) 


which  is  taken  as  single  valued  in  the  complex  plane  apart  from  a  slit  along  the 
negative  real  axis. 

When  the  results  of  Appendix  A  (in  particular  equation  (A-55))  are  applied, 
in  conjunction  with  ( B— 11),  to  equation  (11)  one  obtains  an  expression  for  G(vt) 
which  is  the  asymptotic  expansion  given  in  equation  (3)  plus  an  infinite  number 
of  terms  involving  the  6  function  and  its  derivatives.  Thus  we  obtain 


G(vx) 


w 

■  l 

s=0 


(-)S  ' s  !  N 


(vx) 


s+1 


H(t)  + 


6(t)  +  -j  <5 *  (x )  +  ., 


-  (y  +  log  vx) 


N0  N1 

-f  6(x)  +  -£  6'(t)  +  .., 
v 


\2N1 

— -  <5 (x)  +  ... 

I 

+  . . . 

v 

(14) 


As  pointed  out  in  Appendix  B  such  series  as  those  enclosed  in  the  |  }  in 
(14)  represent  functions  which  are  transcendentally  small  compared  with  the 
negative  power  of  x  .  The  same  is  clearly  true  of  log  x  or  x  m  x  such  a 
series.  Thus  the  term  in  the  square  brackets  in  (14),  which  we  will  call 
G(vx)  ,  satisfies  the  condition 

^JxnG(vx)|  =  0  (n  *  0)  (15) 
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and  so  has  a  zero  asymptotic  expansion  in  terms  of  the  gauge  functions  (vt)”111  . 
Any  function  satisfying  (15)  can  be  interpreted  ( of  Appendix  B)  as  the  sum  of  a 
function  which  is  zero  for  t  greater  than  some  finite  value,  a  number  of  a 
decaying  exponential  expj-  ai^J  (a  >  0,  b  >  0),  and  possibly  a  finite  number  of 
6  functions*.  There  is,  however,  no  reason  why  G(vx)  should  have  non-smooth 
behaviour  at  any  positive  x  . 

One  further  point  to  note  is  that  the  aerodynamic  (or  virtual)  inertia 
matrix  A  ,  is  zero  if  the  flow  is  compressible.  Finally  we  write  the  equation 
of  motion  (1)  in  the  following  alternative  form,  using  (2), 


A  t4  *  v  If +  (E  •  v2)q  =  y2  /  g(vt  *  vv 

dx  y 


dq(T0) 


3  APPROXIMATIONS  TO  THE  EQUATION  OF  MOTION 

The  most  common  approximation  is  to  assume  that  the  function  G  ,  in  the 
expression  for  the  indicial  aerodynamic  matrix  (equation  (2)),  is  zero  and  that 
the  constant  matrices  AQ,  Kq  have  the  values  -B(v),  -  C(v)  respectively 
where  v  is  an  assumed  value  of  the  frequency  parameter.  Thus  instead  of 
equations  (6)  and  (9)  the  approximations 


-  C(v) 


-  B(v) 


are  made.  The  purpose  of  this  approximation  is  to  obtain  a  differential  equation 
which  has  a  solution  which  is  identical,  for  large  time,  with  a  solution  of  the 
original  equation  ((1)  or  (16))  in  the  particular  case  when  the  latter  solution 
ultimately  becomes  sinusoidal  with  frequency  parameter  v  .  This  differential 
equation  has  the  form 


A  §•  +  vb<»)  |a  ♦  (v^C(v)  +  E)q 


*  This  expression,  hereinafter,  refers  to  a  selection  from  6(x)  and  all  the 

«(n)<x)  . 
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A  rather  better  approximation,  due  to  Richardson  is  to  approximate  to  G 
in  equation  (2)  by  a  finite  power  series  multiplied  by  a  decaying  exponential 
term 


G(vt) 


-P0VT 


m- 1  v  ,  ..  r 

£  Kr(VT) 

r=0 


(21) 


Here  m  is  a  positive  integer  and  p^  a  positive  scalar.  In  the  two- 

dimensional  incompressible  case,  for  the  modes  of  heave  and  pitch,  an  established 

9 

good  approximation  to  G  (see  Lomax  )  has  the  temporal  behaviour 


,  -0.09vt  .  „  -0.6vt. 
(e  +  2e  ) 


(22) 


This  suggests  that  a  suitable  choice  for  p^  in  the  approximation  (21)  lies  in 
the  range  0.09  -*  0.6.  With  this  approximation  (equation  (21)) 


K  (p  v)r 

1g(£)  =  V  F--0—  ■ 

v  v  k  (p +  P0V) 


r+1 


(23) 


It  will  be  noticed,  comparing  equation  (3)  or  (14)  with  (21),  or  (11)  with  (23), 
that  this  approximation  to  G  (and  hence  to  the  indicial  aerodynamic  matrix) 
does  not  have  the  right  behaviour  as  t  tends  to  ;  nor,  of  course,  does  its 
transform  have  the  right  behaviour  as  p  tends  to  nought,  or  -  p^v  .  A  suitable 
procedure  for  the  determination  of  the  K^_  coefficients  is  described  in  Ref  3. 

In  particular  from  equations  (9),  (10)  and  (23) 


An  =  -  B 

0  « 


Ko  =  co  ■  c» 


(24) 

(25) 


In  the  applications  of  this  approximation,  in  Refs  3  and  4,  the  relationship  (25) 

was  not  used,  it  being  preferred  to  give  equal  weight  to  the  values  of  the 

oscillatory  aerodynamic  matrices,  at  a  set  of  frequency  parameters,  in  the 

determination  of  the  K  coefficients. 

r 
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4  THE  NATURE  OF  THE  SOLUTION 

4. 1  General 

Taking  the  Laplace  transform  of  equation  (1)  or  (16)  we  obtain  the 
characteristic  equation  which  can  be  written  in  the  two  equivalent  forms 

(A  -  Aj)p2  +  E  -  v2pK(p)  -  0  (26) 

-  Appv  +  (E  -  KjV2)  -  pvG (2)  «  0  .  (27) 

The  number  of  roots  of  this  equation  is  not  necessarily  equal  to  2n  ,  where  n 
is  the  number  of  degrees  of  freedom,  as  it  would  be  if  the  equation  of  motion 
were  a  second  order  differential  equation.  Consider,  for  example,  the  particular 
case  of  two-dimensional  incompressible  flow.  Substituting  equation  (12)  in 
equation  (27)  then  gives 

v2|A(f)  ”  (Aq  +  R)  £  +  “  K,  “  T  +  s)  (  s(£)  “  l)|  J  “  0  •  (28> 

This  function  is  single  valued  in  the  domain  of  the  complex  plane  bounded  by  two 
circles  around  the  origin,  of  very  large  and  very  small  radius,  and  two  lines 
just  above  and  just  below  the  negative  real  axis  -  ie  the  whole  complex  plane 
apart  from  a  slit  along  the  negative  real  axis.  On  the  boundaries  of  this  domain 
2 (2)  has  the  following  values  therefore 

(i)  on  the  large  circle  E  ■  0.5; 

(ii)  on  the  small  circle  E  *  1.0; 

(iii)  on  the  edge  of  the  slit  the  St  (S)  is  always  positive;  and  the  y(S) 
is  always  negative  above  the  slit  and  positive  below  the  slit  except  for  the  two 
ends  where  it  is  zero  (of  Table  1 ) . 

Q 

If  then  we  apply  the  principle  of  the  argument  (Copson  )  to  the  function 
(28)  as  we  traverse  the  boundary  of  this  domain  we  find,  for  the  case  of  one 
degree  of  freedom  (n  ■  1)  that  the  function  can  have  either  two  or  three  zeroes 
within  the  domain.  Consequently  for  n  degrees  of  freedom  the  number  of  zeroes 
can  be  anything  between  2n  and  3n  . 

2 

There  may  in  addition  be  some  negative  real  roots.  Milne  argued  that, 
because  the  imaginary  part  of  the  function  which  is  equated  to  zero  in  the 


characteristic  equation  is  not  single  valued  at  points  on  the  negative  real  axis 
then  the  imaginary  part  cannot  become  zero  at  a  point  on  the  negative  real  axis. 
However,  writing  £r  at  a  point  (£  =  -x)  on  the  negative  real  axis  as 

rr(-x)  =  £r(x)  i  inr(x)  (29) 

then  the  characteristic  equation  (28),  which  has  the  form 

n 

Y  ‘r@  Sr  $  -  0  (30) 

r=0 


is  satisfied,  at  such  a  point,  if 


and 


n 

aQ(-x)  +  ar(-x)^r(x)  =  0 

r=l 


n 

7  ar(-x)nr(x)  =  0 

r=l 


> 


(31) 


both  have  the  same  positive  real  root  (for  x) .  There  is  no  apparent  reason  why 
this  could  not  be  so  and  we  therefore  cannot  rule  out  the  possibility  of  there 
being  a  negative  real  root. 

It  is  not  as  easy  to  determine  limits  on  the  number  of  roots  of  the 
characteristic  equation  for  the  general  case  (not  two-dimensional  incompressible) 
but  the  above  example  makes  it  obvious  that  there  may  be  more  than  2n. 

With  the  use  of  (II)  the  characteristic  equation  (27)  can  be  written 


00  00 

(A-L,,(E)2  -  <vv* *  (4-0  - 1  Ls(*r  4  »•»$  I  «.©' 

'v  '  s=2  s-0 


. (32) 
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At  first  sight  one  may  think  this  equation  has  an  infinite  number  of  roots,  but, 
as  we  have  just  seen  this  is  not  necessarily  so.  One  surmises  that  the  number  of 
roots  will  always  be  finite  and  equal  to  or  a  little  greater  than  2n.  Knowledge 
of  these  roots,  and  the  associated  characteristic  vectors,  will  enable  us  to 
derive  an  expression  for  the  inverse  of  the  characteristic  matrix 

v2M(f)  “  AP2  "  AQpv  +  E  -  K^v2  -  pvG(J)  (33) 

and  hence  the  response  of  the  system  to  an  arbitrary  impulse.  Ue  will  assume  all 
the  roots  of  (32)  are  distinct.  Let  these  be  A^  (i  *  1  -*•  k)  ;  and  a!,  0^  the 
associated  left  hand  and  right  hand  characteristic  vectors.  Then* 


Limit 
p-*A 


l  '  ll  du  I  X^  i 


=  (say) 


(34) 


and  so  M  1  can  be  written,  remembering  (11),  in  the  form 


00  oo 


«"($  -  *  E  E 

i=*l  1  r»0  s=»0 


(35) 


provided  (e  -  K^y2)  is  not  singular.  In  these  equations  the  ou,  0^,  A^  (and 

hence  W.  )  will  either  be  real  or  occur  in  complex  conjugate  pairs;  the  T 

1  ITS 

matrices  will  be  real.  In  particular 


and 


Too  ■  <E  -  V2)"  ♦  7  E  £  <36) 

1=1  1 

T01  "  (E  -  Kov2)"VE  ■  V2)"  * 


(37) 


12 


We  notice  also  that  all  the  Tn  (s  ^  0)  are  zero  if  Nrt  is  zero.  In 

Us  0 

Appendix  A  an  expression,  (A-37),  was  obtained  for  the  inverse  transforms  of 
logmp  .  Taking  this  in  the  form 


6  (x)  =  (-) 

m 


m 


Iso- 


(0)  10gm_VT 


6(t) 


m=l  .  tv 

-ISO- 


(0)  log^-'t 


H(t) 


(38) 


where  calculated  values  of  the  coefficients  w  (0)  are  given  in  Table  3;  it 
follows  that,  for  t  >  0  and  (e  -  non-singular,  the  inverse  transform  of 

M  1  (£)  has  the  form,  using  (B-ll) 


-  *  I  vv  ♦  E  I s-w 


i=i 


*  (say)  X_(t) 


*■  /  Cm  J  (yx) 
5=2  r=0  ^VT' 


(39) 


i  * 


where  the  matrices  are  real.  In  particular  we  note  that 


S02  =  T01 


S12  "  4T02 


S03  "  "  6)T02  "  2TI1 


etc 


(40) 

(41) 

(42) 


and  that  all  the  are  zero  if  (and  hence  -  each  T^s)  is  zero.  The 

general  form  for  the  S  is 
°  rs 


S 

rs 


r+s-1 

E 


u=l 


„(rs)_ 
a  T  • 
u  u-Itr+s-u 


where  the  coefficients 


a (rs) 


can  be  determined  from  (38)  and  (B-ll). 


(43) 
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It  will  be  noticed  that  an  infinite  series  of  terms  in  the  6  function 
and  its  derivatives  has  apparently  been  omitted  from  (39).  The  coefficients  in 
this  series,  which  by  (B-ll)  are  constants,  cannot  easily  be  evaluated  since  they 
are  the  limit  (t  -*■  0)  of  rather  complicated  sums  of  terms  each  of  which 
individually  becomes  infinite  at  the  limit.  We  assume  that  each  coefficient  has 
an  appropriate  finite  value*.  Then,  as  shown  in  Appendix  B,  the  series  omitted 
from  (39)  represents  a  function  which  is  either  zero  for  t  greater  than  some 
finite  value,  or  else  decays,  as  t  -*■  00  ,  as  a  simple  exponential  e  or  more 
rapidly.  There  is  no  reason  why  Xq  should  have  any  non-smooth  behaviour  at 
finite  x  ,  and  we  have  already  included  in  (39)  all  the  possible  simple 
exponential  terms.  Consequently  the  function  omitted,  if  it  is  not  a  null 
function,  must  be  a  function  which  decays  more  rapidly,  with  increasing  x  , 
than  any  of  the  terms  we  have  given  in  the  above  expression  for  Xq 
(equation  (39)). 

An  arbitrary  instantaneous  excitation  at  r  «*  0+  can  be  taken  to  be 
6(T)fQ  +  6(1)(x)f]  +  <5(2)(x)f2  +  ... 

where  fQ,  fj  ,  ...  are  arbitrary  real  constant  column  matrices;  and  so  the 
free  motion  of  the  system  (ie  the  motion  after  the  disturbance)  is  given  by 


where 


00 


(44) 


(45). 


and  Xq  is  given  by  equation  (39). 

Milne  has  derived  the  same  form,  for  the  fundamental  solution  Xq  ,  in 
Ref  2  by  a  different  method.  The  infinite  of  arbitrary  constants  (the  elements 


*  Thus  the  coefficient  of  6(t)  contains  —  *  a  power  series  in  (log  x/x),  and 
in  the  simplest  case,  for  example,  we  can  make  the  identification 


I 


v*  /,„  ,\'l 


/ 


\ 


Sf t 
x-*0 


]_ 

x 


S't 

x-K) 


-  0  . 
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of  Che  matrices  fj  )  arises  from  Che  fact  that  Che  free  motion  subsequent  Co  a 

given  instant  depends  not  only  on  the  velocity  and  displacements  at  that  instant 

(as  with  an  instantaneous  dynamical  system)  but  on  the  whole  history  of  the 

motion  previous  to  that  instant.  For  large  time  the  ~  terms  in  the  above 

solution  will  be  dominant*;  for  intermediate  times  the  exponential  terms  need 

also  to  be  taken  into  account,  and  indeed  Milne  has  suggested,  with  some 

supporting  evidence,  that  in  practice  the  leading  S  coefficients  (eg  (39)) 

r  s 

are  small  compared  with  the  and  so  there  will  be  a  large  range  of  time  for 

which  the  motion  is  predominantly  exponential.  For  very  small  time  the  form  of 
solution  that  we  have  obtained  is  of  course  inadequate. 

4.2  Particular  cases 

At  this  point  we  will  note  certain  features  of  the  fundamental  solution 

which  appears  in  particular  cases.  Firstly  if  the  flow  is  three-dimensional  the 

matrix  N_  ,  and  hence  each  matrix  S  „  ,  is  zero  and  so  the  fundamental  solution 
0  r2 

has  the  form 


l  vv  -3££  u* 


s=3  r=0 


in  which  the  dominant  term  for  large  time  is  of  Of-L)  .  Secondly  if  the  matrix 
(e  -  K^v  )  is  singular,  that  is  if  the  characteristic  equation  (32)  has  a  zero 
root  Aj  *  0  ,  then  the  dominant  non-exponential  term,  apart  from  the  constant 
term,  is  of  0(^-)  (of  Ref  6)  and  so  the  fundamental  solution  has  the  form 


0  “  V(W1  +  Yj  ^  1  )  +  V  Z  Yj  l0g2(VT) 


U  _(v) 


s= 1  r«0 


This  can  be  simply  demonstrated  for  the  one  degree  of  freedom  system  by  expanding 
in  the  form  of  equation  (35),  obtaining  its  inverse  transform  using  the 

v  -I 

results  of  Appendix  A  and  hence  obtaining  the  inverse  transform  of  M  .  The 
two-dimensional  case,  when  the  characteristic  equation  has  a  zero  root,  is  not  as 
simple. 


*  Assuming  all  the  have  negative  real  parts. 
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We  take  as  our  independent  variable  P  »  which  is  a  monotone 

function  of  p  ,  and  note  that 


Limit 

P^A. 


X. 


P  l°g(£)  A. 


v  -  p 


v  -  X . 
i 


M-1® 


6. a! 

l  l 


at/ — dM(u) — \  „ 


>  (48) 


W*  (say) 


where  a!,  3.  are  the  left  hand  and  right  hand  characteristic  vectors  associated 
i  i 

with  the  root  p  =  X^  .  If  i  =  1  denotes  the  zero  root,  ie  X^  »*  0  ,  then  for 
i  ^  1  the  W*  are  related  to  the  W^  (equation  (34))  viz 


W*  «  Limi  t 
1  P^Xi 


vW 


j.  /P  l08  (y)  Ai 


P  ~  Xi  \  v  -  p  v  -  X. 


(49) 


Consequently  the  inverse  of  M  can  be  expressed  from*  (48)  and  (49)  in  the  form 


«•'  © 


C>  -  9 


CO  00 

\  ^  \  ’  rTi\ 


T*VL  T^TJ  "2Ll  $"*  ■  (50, 

1  r=0  s=0 


^  log  —  .  „ 

v  v  1*2 


This  gives  a  fundamental  solution  (the  inverse  transform  of  M  *  )  of  the  form 
(see  Appendix  A) 


-  vW* 


y  (-,« _ j£1__  ,  x  y 

L  (x  +  iog  vr)U+1  L 


(-) 


U  (u  +  »wu(x) 


(X  +  log  vt) 


u+2 


"iv1  * y3  Yj  l  io®r 

i*2  s*2  r*0 


(vt) 


S*  ,  V 
rs(v) 


(vt) 


r+s 


for  t  >  0 


y  (so 


*  We  are  still  assuming  all  the  X.  are  distinct.  In  addition  we  assume  that 
no  ■  v  ,  otherwise  equation  (50)  will  require  some  modification. 
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for  Che  present  case  when  the  flow  is  two-dimensional  and  the  characteristic 
equation  has  a  zero  root.  The  wy  coefficients  are  given  in  Table  3,  having 
been  determined  from  equations  (A-34)  and  (A-35).  The  value  of  x  is  arbitrary, 
but,  as  explained  in  Appendix  A,  for  any  particular  value  of  vr  ,  it  has  to  be 
suitably  chosen  so  that  a  sufficiently  accurate  'sum'  can  be  obtained  from  the 
asymptotic  series.  Thus,  for  example,  for  vt  =  1  ,  with  x  =  -3y  >  where  y  is 
Euler's  constant,  the  first  series  can  be  evaluated  to  about  three  significant 
figures;  but  with  x  *  3y  not  a  single  significant  figure  is  obtained  ( of 
Appendix  A).  For  r  large,  from  equations  (A-45)  and  (A-53)  or  otherwise 


oo 

E 

u=0 


(->l 


VU(X> 


(X  +  log  vt) 


u+1 


1 


log  VT 


(52) 


Similarly 


OO 

E 

u=0 


( ~) 


(u  ♦  \  >wu  (X) 


(X  +  log  vt) 


u+2 


,  2 
log  VT 


for  t  large 


(53) 


Consequently  in  this  case  (two-dimensional,  characteristic  equation  with  zero 
root)  the  dominant  term  for  large  time  is  Q(f0'g  ~T)  provided  all  the 
A^  (i  =  2  -*■  k)  have  negative  real  parts.  The  T*g  matrices  appearing  in  (50) 
are  given  by  rather  complicated  expressions.  In  particular 


T* 

00 


v2NqW*  +  (E  -  Kv2) 


K. 

E 

i»2 


E? 

i-2 


-1 


:  (wi +  Wi +  wiNowT) 


W.  (E  -  K  v2) W. 

+  V  -ii _ 2—1 -1 

/.  A. A. 

f—i  l  j 

J*2 


> 


(54) 


The  relationship  between  the  S*  matrices  and  the  T*  matrices  is  the  same  as 

rs  rs 

that  between  the  Srg  and  the  Trg  (equations  (40)  to  (43)).  It  is  perhaps  a 
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little  surprising  to  find  ( of  equation  (45)) 


<St  (X.)  =  0  for  all  j  (55) 

T-*oo  J 


in  a  case  where  the  characteristic  equation  has  a  zero  root.  One  would  expect 
that  there  would  be  excitation,  which  ultimately  died  away,  which  would  produce  a 
response 

q  *  8H(t)  (56) 


That  is  there  would  be  a  means  of  producing  a  'free'  motion  which  was  a  constant 
displacement.  However  if  we  surmise  the  displacement  (56)  then  from  equation  (14) 
the  necessary  excitation  is 


Since 


k 


A5"(t)  - 


V6 


(t)  - 


v2G(vt) 


(« -  v2)) 


BH(t) 


£ft  G(vt)  =  0 

X"H» 


(57) 


(58) 


then,  as  r  tends  to  ®  ,  this  force  tends  to 


(E  -  Kcv2)B  . 


(59) 


In  the  case  under  consideration  (e  -  K^v2)  is  singular  and  so  the  excitation 
force  (57)  will  ultimately  tend  to  zero  if  8  in  the  right  had  eigenvector  8( 
of  (e  -  K^v2)  .  A  'free'  motion  of  constant  displacement  is  therefore  possible. 
Such  motion  cannot  however  be  achieved  by  an  instantaneous  excitation  for  we  have 
from  (14) 


G(vt)8, 


uu 

y  irOf : 

4*  (VT) 


S+l  1 


N  8,H(t) 


s+l  s^l 

_  I  1»T  1 

s=0 

+  an  impulsive  term 


(60) 


For  this  reason  we  have  used  inverted  commas  in  designating  such  motion  as  'free' 
motion.  In  the  three-dimensional  case,  when  the  characteristic  equation  has  a 
zero  root,  it  follows  from  (47)  that  an  impulse  6(T)f^  will  produce  a  response 
which  is  ultimately  a  constant  displacement  but  here  again  a  response  of  the 
form  (56)  can  only  be  achieved  by  an  excitation  which  is  not  entirely  impulsive. 
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4.3 


Discussion 


To  recapitulate  we  see  that  we  have  obtained  four  forms  for  the  funda¬ 
mental  solution  Xq  .  If  the  matrix  (e  -  K  v^}  is  not  singular,  that  is  if  the 
characteristic  equation  does  not  have  a  zero  root,  is  given  by  equation  (39) 

or  (46)  according  as  the  flow  is  two-dimensional  or  three-dimensional.  When 
(e  -  K^v2)  is  singular,  however,  we  have  the  modified  forms  (51)  and  (47)  for 
two-  and  three-dimensional  flow  respectively.  In  every  case  we  have  assumed  that 
the  roots  of  the  characteristic  equation  are  distinct.  Cases  of  multiple  roots 
will  produce  some  further  modification  to  the  fundamental  solution.  The  free 

motion  of  the  system,  for  t  >  0  ,  is  then  given  by  equation  (44).  It  will  be 
.  .  .  X^t 

noticed  that  m  each  case  this  solution  involves  constituents  B^e  ,  though 

there  are  other  constituents  as  well.  One  expects  this,  but  it  should  be  noted 

that  the  number  of  these  exponential  terms  may  well  be  more  than  2n  ,  where  n 

is  the  number  of  degrees  of  freedom,  and  also  that  the  right  hand  characteristic 

vectors  B^  have  no  orthogonal  relationship  with  the  left  hand  characteristic 

vector  a!  .  If  g(^)  were  zero,  as  is  assumed  in  the  first  approximation 

considered  in  section  3  (equation  (19)),  then  {of  equation  (27))  these  vectors 

would  satisfy  the  biorthogonal  relationships 


[°i  Vi] 
[ai  Vi] 


K  v  -  E 
a 


K  v  ~E 
a 


K  v  -  E 

o 


V 


B. 

J 


X .  B . 
J  J 


X .  8  • 
J  J 


=  0 


(i  =£  j)  (61) 


=  0 


(i  #  j)  (62> 


but  there  is  nothing  corresponding  when  no  approximation  is  made. 

5  CONCLUDING  REMARKS 

It  follows  from  the  forms  of  fundamental  solution  derived  in  the  previous 
section  (equations  (39),  (46),  (47)  and  (51))  that  an  aeroelastic  system  will  be 
stable,  in  the  sense  that  the  response  to  any  impulsive  excitation  ultimately 
dies  away  to  zero,  for  small  perturbations  in  subsonic  flow  provided  all  the 
roots  of  the  characteristic  equation  ((26),  (27))  have  negative  real  parts*. 

*  We  have  seen  that  in  two-dimensional  flow  we  can  still  have  stability  in  this 
sense  when  the  characteristic  equation  has  a  zero  root. 
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Thus  the  traditional  UK  method  of  flutter  investigation  with  lined-up  frequency 
3 

parameter  is  completely  adequate  for  determining  such  stability;  for  it  deter¬ 
mines  any  roots  with  zero  real  part  correctly  and  so  by  a  survey  of  the  speed 

,  4 

range  finds  whether  any  root  becomes  unstable.  There  is  also  some  evidence  that 
it  determines  complex  roots  with  a  fair  degree  of  accuracy  even  when  the  real 
part  is  negative  and  relatively  large  (compared  with  the  imaginary  part) . 

Often  one  wishes  to  know  not  only  the  critical  flutter  speed  of  an  air¬ 
craft,  if  there  is  one  in  the  flight  domain,  but  also  how  quickly  the  effect  of  a 

£ 

disturbance  will  become  insignificant.  If  Milne's  supposition  that,  in 

equation  (39)  for  example,  the  leading  S  coefficients  are  usually  small 

IT  S 

compared  with  the  coefficients,  ,  of  the  exponential  terns,  then  the 

traditional  UK  approach  should  also  provide  this  information  reasonably  well.  The 
characteristic  value  which  has  the  least  negative  real  part  given  the 

required  information. 

However  when  the  non-exponential  terms  in  the  fundamental  solution  are  of 
significant  size  at  intermediate  times  one  would  have  to  consider  also  at  least 
the  dominant  one  of  such  terms,  which  in  the  two-dimensional  case  is 
(see  equations  (37),  (39)  and  (40)) 

4  (E  -  Kav2)_1VE  -  V2)'1  (63) 

T 

and  in  the  three-dimensional  case*  (see  equations  (42)  and  (46))  is 

-  4  (E  -  K/fX  (E  -  Kov2)_1  .  (64) 

T 

If  one  determined  the  appropriate  expansion  of  the  oscillatory  aerodynamic 
matrix  C(v)  +  ivB(v)  for  small  v  then  ( of  equations  (8)  and  (11))  this  would 
provide  the  required  matrix  or  Nj  . 

In  this  note  we  have  only  considered  the  free  behaviour  of  our  aeroelastic 
system.  The  fundamental  solution  Xq(t)  could  in  principle  be  also  used  to 
obtain  the  response  to  any  excitation  by  evaluating  the  faltung  integral  of  Xq 


*  Since  N, 


is  zero  so  is 
T 

*11 


02 


and  it  is  found  then  that 


(E  -  K^f'NjfE  -  Kov2)-' 


in  a  form  convenient  for  use  at  small  t  which  we  have  not  obtained.  It  would 
probably  therefore  be  advisable  to  determine  first  the  Laplace  transform  of  the 
response  and  then  invert  that. 
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CO 

r.  (n)  /  x  f  - u  ,  n  2-1 
r  (z)  =  I  e  log  u  u 


du 


(A-2) 


and,  putting  z  ■  1  ,  we  have 


00 

r^(l)  ■  /  e  u  lognu  du  . 


(A-3) 


The  Laplace  transform  of  logmt  ,  for  m  >  0  ,  is  therefore  given  by 


00  00 
J  e  pt  logmt  dt  m  p  J  e  UU°8  u  “  log  P 


du 


(A-4) 


m 


r(r>(D  iog"-rP 


13,16 


Wong  ’  has  shown,  inter  alia,  that  a  similar  formula  applies  for  negative 
powers  of  log  t  .  He  considers  the  transforms  of  functions  f(t)  which  are 
locally  integrable  and  satisfy 


f(t)  ~  log  nt 


as  t  -*■  00  (A-5) 


where  n  is  positive,  and  so  he  proves  in  particular  that  the  Laplace  transform 
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of  log  nt  has  the  asymptotic  expansion* 


e  pt  log  nt  dt 


^  £  0  -'-Ssi 


as  p  -*■  0  .  (A-6) 


(This  is  an  expansion  using  the  sequence  of  gauge  functions 

Jl-L—  _ 1  \ 

.  n  ’  .  n+1  ’  ***  j  • 

p  log  p  p  log  p  / 

The  values  of  the  derivatives  of  T(z)  at  z  =  1  can  be  obtained  from 
the  well  known  series^  for  -r-  {log  r(z)|  ,  ie 


y  (z)  =  r  ^ 

Vz;  r(z) 


-  -y-l.f 

z  yn  n  +  zj 


(A-7) 


where  y  (=  0.577215....)  in  Euler's  constant, 
^(z)  and  its  derivatives 


oo 

T  (z)  =  {log  r(z)}  =  r(r)(-)r  y  - ! — - 

dz  f-A  (n  +  Z) 


(A-8) 


are  the  polygamma  functions.  A  simple  general  expression  for  the  rw(z)  cannot 
be  written  down.  However,  differentiating  r^(z)  |=  1'j(z)r(z)|  (n  -  1)  times, 
using  Leibniz's  theorem,  we  obtain  the  recurrence  relationship 


r(n)(z) 


£  CD 


(z)  , 


(A-9) 


and  so,  in  particular 


r(n)(l) 


Zr(n)  „ (n-r ) 

r(n  -  r  +  1)  prJ 


(A- 10) 


*  The  meaning  of  this  is  that  the  sum  of  a  finite  number  of  terms  (N)  differs 


from  the  required  function  by  —  of - j  for  every  integer  N  £  0  . 

p  \logn  p/ 
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where 

and 


*,0) 

1  =  TTiT 


»r(D 


=  -  Y 


r(r) 


(A- 1 1 ) 


l,  7 


(r  *  2)  .  (A- 12) 


A  more  general  relationship  than  (A- 10),  involving  the  operator  (D  +  x)  ,  where 

D  =  ,  will  however  be  useful  to  us.  To  obtain  this  we  note  that  Leibniz's 

dz 

theorem  can  be  generalised  to  the  form 


n 

(D  +  a  +  B)n(uv)  =  ^  ^  J  (D  +  a)r (u) (D  +  B)n  r(v) 


(A- 13) 


of  which  a  particular  case  is 


/  n- 1  \ 

(D  +  x)n_,{'*',  (z)r(z)|  =  V  (  j '?r+1(z)(D  +  x)D'1"r{r(z)} 

r=0  ^  r  ' 


(A- 14) 


where  x  is  arbitrary. 
Thus  we  have 


n  .  j . 

(D  +  x)n{r(z)f  -  y  (  J  (z) (D  +  X)n“r{r(z)[  +  X (D  +  X)n_1 |r(z)f  (A-I5)- 
r=l  r  1 


and  putting  z  «  I  and  writing 


En(x)  -  [<D  ♦  x)n|r(z)|]  -  $  (  ‘  ) 

s«0 


xsr(n-s)(i) 


(A- 16) 


this  gives 
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n 

5„<X>  *  £  r<n  -V*  l)  °r«„-r<»  <”  »  '> 


r=l 
1 


(n  *  0) 


V  (A- 17) 


where 


a,  =  Pj+X  =  X-y 


a  =  p 
r  r 


(A- 18) 

(r  i  2)  .  (A- 19) 


We  now  require  the  inverse  transforms  of  logmp  which  we  will  denote  by 


0  (t)  ,  where 
in 


oo 

J  e  pt0m(t)dt  =  logmp  . 


(A-20) 


From  equation  (A-4) ,  for  positive  m  ,  using  (A— 16) 


OO 

h'pt 


(X  +  log  t)mdt  * 


m-s 


•  £  "*0“*^  £  r™">> — (7)  >  «■”> 


s=0 


r=0 


■  7  1  (-)SC)?»-s<x>  •*' 

s*0  X  s  7 


Let  X  be  the  infinite  lower  triangular  matrix  whose  non-zero  elements 


are 


.  n«-  • 

ij  \j-i )  i-j 


(X) 


(i  *  j)  .  (A-22) 


Thus,  in  particular,  of  equation  (A-17),  X  has  unit  elements  on  the  principal 
diagonal*.  Its  inverse  will  then  have  the  form 


*  If  we  chose  X  ■  y  it  would  (see  (A- 18))  also  have  zero  elements  on  the  first 
subdiagonal. 
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X_1  =  I  +  S  (A-23) 

where  S  is  a  lower  triangular  matrix  with  zero  elements  on  the  principal 

diagonal*.  We  can  then  solve  the  infinite  set  of  equations  (A-21)  (m  =  0  -+  00 ) 

in  terms  of  the  elements  s  of  S  .  Thus 

rs 


co  g  00 

J  e  pt(X  +  log  t)Sdt  +  ^  sg+1  v(X)  J  e  pt(x  +  log  t)V  *dt  . 

0  v—  I  0 

for  s  ^  0  . (A-24) 

It  follows  that  the  inverse  transform  of  logmp  ,  for  m  positive,  is 

I*"  m 

(X  +  log  t)m  +YjU+  log  t)V”lsm+l,v(X)  6(t) 

m 

+  m(X  +  log  t)™  1  +  ^  (v  -  1)(X  +  log  t)V  2sm+,  v(X) 

v=2 


At  first  sight  there  may  seem  to  be  an  arbitrariness  about  this  expression  since 

X  is  arbitrary.  However  it  will  be  apparent,  when  we  obtain  expressions  for  the 

suv(x)  ,  that  this  is  not  so.  Indeed  the  main  purpose  of  introducing  an 

arbitrary  x  ,  rather  than  taking  it  *  y  or  0,  is  in  the  subsequent  derivation  of 

0  (t)  when  m  is  negative, 
m 

From  the  definition  of  S  ,  equation  (A-23),  it  follows  that  its  elements 
along  a  diagonal  are  related 


(A-25) 


*  If  we  chose  X  ■  Y  it  would  (see  (A- 18))  also  have  zero  elements  on  the  first 
subdiagonal. 
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The  first  column  of  (I  +  S)  is  thus,  using  (A-23) 


t“0  ”l  “2  •••! 


where 


W  as  1 
0 


(A-27) 


(A-28) 


n 

“n00  ’  ‘ 


(A- 29) 


Equations  (A-29)  and  (A— 1 7 )  are  however  probably  not  the  most  convenient  way  to 
determine  the  w^(X)  ■  It  is  obvious  that  there  is  also  the  reciprocal 
relationship 


5"<X)  ‘  '  I  Q  Vt(X)”r<X) 


(nil)  .  (A-30) 


Substituting  from  (A-16)  in  this  equation  and  reversing  the  order  of  summation  on 
the  right  hand  side  then  gives 


(1)  (A-3 1 ) 


(1)  =  0  .  (A-32) 


This  equation  is  true  when  n  is  any  positive  integer  and  so  it  is  easily  seen 
by  induction  that  we  must  have 


r(s)(i) 


•  1 0  •-«>'“ 


(i)  =  (-  xV 


(A-33) 


which  can  be  rewritten,  bearing  in  mind  (A-28),  as 
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“n<X)  .  (-  X)"  r  I 

f  Vr(X)  'I 

\(  r(r)(i)  \ 

r(n  +  1)  r(n  +  1)  ' 

i^r(n  -  r  +  1)J 

'  \ r(r  +  1) ) 

(A-34) 


Equation  (A- JO)  gives 


r(n)(i)  r^n~r) (i) 

r(n  +  1)  “  n/,Prr(n-r+l) 

r=l 


(A-35) 


and  so  these  two  equations,  (A-34)  and  (A-35),  provide  a  suitable  means  of 
successively  evaluating  the  (wn(X)/r(n  +  I))  from  the  values  of  the  Pf 
(equations  (A— 11)  and  (A- 12)).  We  have 


P 


r 


r(r)(i) 
r(r  +  1) 


r  -*•  ”  (A-36) 


Values  of  these  two  functions  are  given  in  Table  2.  The  were  evaluated,  to 

rather  more  significant  figures  than  shown  in  the  table,  using  Euler's  summation 
14 

formula  .  Values  of  the  wr(X)  ,  calculated  from  (A-34),  are  given  in  Table  3 
for  several  values  of  X  . 

Equation  (A-25) ,  for  the  inverse  transform  of  logmp  ,  for  m  positive, 
can,  using  (A- 26)  and  (A-28) ,  be  rewritten  as 


0  (t) 
m 


(-)m 


Wy(X) (X  +  log  t)“~V 


6(t) 


+  m 


wv(X)(X  +  log  t)°"v_1 


H(t) 

t 


(A-37) 


Now,  substituting  for  £n(X)  from  equation  (A- 16)  in  (A-30)  gives,  after  a 
little  manipulation 


i 
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while  from  (A-33)  we  have 


Comparing  these  two  equations  we  immediately  see  that  we  must  have 


XS'rwr(X)  *  wg(0)  . 

It  follows  therefore  that  there  is  no  arbitrariness  in  equation  (A-37) . 

Let  us  now  consider  the  inverse  transform  of  logmp  ,  when  m  is 
negative.  From  (A-6) ,  using  (A- 16)  and  (A-22) ,  since  we  can  add  asymptotic 
series  as  if  they  were  convergent, 


00 

I 


rpt 


(X  +  log  t) 


n 


dt 


(->n  y  rn\  (-)Vrj(o 

P  r=0  '  r '  ^log  P  “  x>n+r 

nJiLLly  r  w><,> 

k  i°*n+8p  r4*w 

(->-  y  <-)S  (?)  {s(x) 


as  p  0 


±L> 

P 


,  n+s 
log  p 


(n-lK-n 

oo 

/  \n  V-'  X 
(-)  \  _ rn 

P  Lj  ioor 


(X) 


r«n 


log  p 


> 


(A-38) 


(A-39) 


(A-40) 


(A-4 1 ) 
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This  result  can  be  confirmed  by  the  use  of  Wang's'^  general  formula  for  the 
Laplace  transforms  of  function  of  the  form  (A-5) .  Consequently 


OO 

<->"  f 


(X  +  log  t) 


dt  + 


OO  On 

E.  ^ 


rpt 


v=n+1 


Vn  J  (X  +  log  t)V 


dt 


y is.,  y  .  yjk 

/  i  i  r  /j  vn  /  j  .  r 
log  p  i 

r=n  &  K  v*n+l 


r=v 


log  p 


as  p  ■+  0 


> 


oo  oo  Y* 

yj^.  y  _!_  y  , 

logrp  ,  log  p  L—> . 

r=n  °  v  r=n+ 1  v 


s 

rv  vn 


v=n+ 1 


i  n 
P  lc>g  P 


J 

(A-42) 


since,  from  (A-23) 


r 

z 

v=n+l 


rv  vn 


-  x 


rn 


(r  >  n) 


(A-43) 


and  xRn  is  unity.  Since  the  asymptotic  expansion,  which  in  this  case  termin¬ 
ates  after  just  one  term,  is  in  terms  of  the  gauge  functions  ( - - — ) 

p  logn  rp 

(r  ■  0  -*■  «>),  it  follows  that  the  difference  between  the  left  and  right  hand  sides 

of  (A-42)  is  transcendentally  small,  as  p  ■+■  0,  compared  with  these  gauge 

functions.  That  is  this  difference  is  of - — r — 1  for  any  r  however  large; 

p  logn+rp 

or,  putting  it  another  way,  it  is  a  function  which  has  a  zero  asymptotic 
expansion  in  terms  of  these  gauge  functions.  It  follows,  therefore,  with 
similar  significance,  that,  using  (A-26)  and  (A-28) ,  for  n  positive 
{of  equation  (6-12)),  we  have  the  following  asymptotic  expansion,  which  also 
always  terminates  after  one  term. 
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I  du  (X,t) 

SS  P  (X,0)<S(t)  +  — ^rr -  H(t)  ~  — 

In  dt  ,  n 

I  log  P 


as  p  -*•  0  (A-44) 


where 


Z/n+v-1 \  w 

(-)V[  ) - - 

v=0  V  v  J  (X  +  lc 


w  (X) 


(A-45) 


log  t) 


Examples  of  terms  which  are  transcendentally  small  compared  with  the  gauge 
functions  log  ^n+r^p  ,  as  p  +  0  ,  are 


P  (s  >  1) 


(p  +  a) 


(s)  ~3.t. 

which  have  inverse  transforms  6  (t)  ,  <$(t)  -  ae  respectively.  These 

latter  functions  are  transcendentally  small  compared  with  the  gauge  functions 
(X  +  log  t)  ,  as  t  4  ”,  If,  as  seems  reasonable,  we  can  assume  that  this  is 
also  true  for  the  inverse  transforms  of  any  functions  which  are  transcendentally 
small  compared  with  the  log  ^n+r^p  ,  as  p  4  0,  then  we  can  write 
(of  equation  (A-20)) 


6  (t)  ~  y  (X,0)4(t)  + 

- Tl 


dy  (X , t) 


H(t)  as  t  4  »  (A-46) 


or  equivalently 


dy 

e-n(t>  ~  S f 


as  t  4  ”  .  (A-47) 


Thus,  substituting  from  (A-45),  we  have  the  following  asymptotic  expansion  which 
can  be  used  to  evaluate*  the  inverse  transform  of  log  np  ,  when  n  is  positive. 


9_n(t)  ~  (-)* 


ft 

h  v  v )  <x  * 


w  (X) 

V 

log  t) 


n+v+1 


t  -*■  “  .  (A-48) 


*  With  an  asymptotic  expansion,  if  the  series  is  divergent,  a  point  will  be 
reached  where  the  addition  of  further  terms  increases  the  error.  One  must 
therefore  stop  at  this  point. 
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It  should  be  noted  that  there  will  be  a  limit  on  the  arguments  for  which 
this  equation  is  of  use.  The  error  might  contain  some  decaying  exponential  terms 
and  so,  even  though  the  series  on  the  right  of  (A-48)  will  converge  for 
sufficiently  small  t  ,  the  value  it  gives  for  6_^  may  be  appreciably  in  error. 
However  there  is  no  reason  why  the  expansion  should  not  be  used  for  negative 
values  of  (X  +  log  t)  and  indeed  it  may  prove  more  convenient  to  take  a  value 
of  X  so  that  this  is  so. 

As  an  illustration  of  this  let  us  consider  the  associated  function 


00 

^  (t)  =  -  [  e 

n  J  -n 


(x)dT 


(A-49) 


v  (t)  ~  p  (X,t) 

n  n 


(A-50) 


Taking,  for  example,  a  positive  value  of  X  and  trying  to  use  this  asymptotic 
expansion  to  evaluate  Vj(l)  it  is  found  to  be  of  no  practical  use.  Thus  with 
X  =  3y  successive  sums  are: 


Number  of  terms 


-0,58  -0.96  -0.90  -0.17  1.17  2.33  1.30  7.65. 


and  for  larger  X  things  get  worse  (cf  the  behaviour  of  the  w^(X)  as  shown  in 
Table  3).  However  if  we  take  X  negative  there  is  no  difficulty.  Thus  with 
X  =  -  3y  the  successive  sums  for  Vj(l)  are: 


Number  of 
terms 


1  2  3  4  5  6  7  8 


20  21 


0.58  1.35  2.06  2.43  2.40  2.23  2.19  2.28  _  2.29  2.22 


Evidence  that  the  sign  of  X  does  not  make  any  change,  such  as  the 

addition  of  a  constant,  to  the  set  of  functions  for  which  p  (X,t)  is  an 

n 

asymptotic  expansion,  is  given,  for  example,  by  consideration  of  the  evaluation  of 
Vj  (t)  for  t  large.  Let  us  take 


X  ■  m  log  t 


(A-51) 
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where  m  is  any  number,  positive  or  negative,  othev'  than  -1  such  that  |X|  will 
also  be  large.  Then,  when  that  is  so  (of  equation  (A-34))  we  have 


wy(X)  **  (-  X)  . 


(A-52) 


It  is  then  easily  seen  that  the  sum,  of  N  terms  of  the  series  for  Pj  (X,t)  , 
equation  (A-45),  and  hence  of  the  asymptotic  expansion  for  v^(t)  ,  is 


vj(t)  ~  Vx’0  *  “  Tog 


(A-53) 


Thus  we  ultimately  get  the  same  value  whether  m  is  positive  or  negative  (or 
indeed  zero,  as  (A-45)  and  Table  3  show)  provided  —  -yy|  <  1  . 


<  1  . 


We  have  thus  established  the  validity  of  our  expression  (A-48)  for  the 
inverse  transform  of  lognp  when  n  is  a  negative  integer.  It  should  perhaps 
be  noted  that  we  have  there  an  infinity  of  asymptotic  expansion,  for  the  para¬ 
meter  X  is  arbitrary.  For  n  positive  we  had  previously  obtained  the 
expression  (A-37)  -  again  containing  an  arbitrary  parameter  X  .  To  complete  the 
picture  we  have  the  well  known  fact  that 


eQ(t)  =  6(t) 


(A-54) 


it  being  the  function  whose  Laplace  transform  is  1 . 

A  particular  case  which  we  require  in  the  main  part  of  the  Memorandum  is 
that  when  n  =  1  which  from  (A-37)  is 


e  j  (t) 


-  (y  +  log  t)S(t)  - 


(A- 55) 
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THE  RIGHT  HAND  DIRAC  DELTA  FUNCTION  AND 
ITS  USE  IN  LAPLACE  TRANSFORM  THEORY 


Since  the  Laplace  transform  integral  has  a  lower  limit  of  zero  one  has  to 
use  what  is  called  a  right  hand  Dirac  delta  function  <$(t)  which  is  defined  to 
be  such  that 


6(n)(0)  =  0 
5(n)(»)  =  o 


n  =  0  -*■  °°  (B-l ) 


t 

l  6 (t) f (x)di  =  f (O)H(t) 
0 


(B-2) 


where 


H(t)  =  0 


t  «  0 


*  1 


t  >  0 


(B-3) 


and  f(t)  is  a  fairly  good  function 
It  follows  that  we  can  write 


5<t)  = 


dH(t) 

dt 


(B-4) 


for 


t 

/ 


f(t>  dr 


'  [f<T,H(t>]o  ’  / 


df  (t) 
dt 


H(x)dx 


H(t) f  (t) 


■i 


df  (x) 
dx 


dx 


(B-5) 


H(t)f (0)  +  { H(t)  -  l{{f(t)  -  f(0)} 


H(t)f  (0) . 
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Differentiating  (B-2)  we  see  also  that  6(t)  is  such  that 


6(t)f(t)  =  6(t)f (0) 


(B-6) 


For  the  derivatives  of  6  we  have,  integrating  by  parts  and  using  (B-l), 


6^  (r)f  (x)dx 


n-1 

'  I 

r=0 


(-)rf(r)(t)6(n_1"r)(t)  +  (-)nf (n)(0)H(t)  .  (B-7) 


Also,  from  (B-6),  by  successive  differentiation* 


.(n) 


(t)f(t)  =  ^  M  (-)r6(n  r)(t)f 

r^O 


(r) 


(0) 


(B-8) 


and  so  (B-7)  can  be  rewritten  as** 


t 

/ 


6(n)(x)f(x)dx  =  ( 


n-1  .  . 

l  O 


)r6(r)(t)f(n“1_r)(0)  +  (-)nf(n)(0)H(t)  . 


(B-9) 


It  follows  from  (B-2)  and  (B-9),  putting  f(x)  -  e  pT  ,  t  =  “  ,  and  using 
(B-l),  that  the  Laplace  transforms  of  the  6  functions  are 


^(6(n)(t))  =  J  e'ptS(n)(t)dt  =  pn  n  =  0  +  »  .  (B-10) 

0 

Consequently,  from  the  well  known  theorem  for  the  inverse  transform  of  a  product, 
we  find  that,  using  (B-9)  and  (B-7),  for  n  >,  1 


* 


** 


Incidentally  putting  t  =  0  in  (B-8),  we  notice  that  this  equation  would  not 
be  satisfied  if  the  first  of  the  conditions  (B-l)  in  the  definition  of  6  were 
relaxed. 


We  have  used  the  fact  that 


n-  1-r 


l  O  •  O 
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L 

</  1  {pn£”{f  (t))|  =  J 


f  (t  -  T )  6  (x)dX 


■\ 


to 


6(r)f(n-1-r)(tj  +  fn(tjH^cj 


> 


n- 1 

V  £<r)(0)6<n-'-r)(t)  *  f <n) (t)H(t)  . 
ki  J 


(B-l  1 ) 


The  particular  case  n  =  1  is  the  well  known  formula 

^_1|p^(f(t))^  =  f(°(t)H(t)  +  f  (0)  <5  (t)  (B-l  2) 

which  can  be  written,  using  (B-10),  as 

P^(f(t))  -  f(0)  -  ^|f(0(t)H(t)l.  (B-l 3) 

In  practice  one  may  well  obtain  as  an  inverse  transform  a  function  in  the 
form  of  an  infinite  series  of  terms,  containing  the  5  function  and  its 
derivatives .  At  the  two  limits  t  =  0,  “  such  a  function  is  by  definition  zero, 
but  how  do  we  interpret  it  elsewhere?  Consider  for  example  the  function 


X(t) 


E 

s»0 


(-)s6(s)(t) 


s+l 

a 


[&(&)  >  0)  (B-l 4) 


which  has  the  same  Laplace  transform  as  e  for  p  >  0.  Lerch's  theorem  states 
that  if  two  functions  have  the  same  Laplace  transform  then  their  difference  is  a 
null  function*.  Thus  for  this  example  we  require  therefore  that 
t 

J  |e"aT  -  X(x)}dx  ■  7  {»  -  e‘at  +  *(t)  "  «(t)|  (B-15) 

0 


*  A  null  function  is  one  whose 


/ 


is  zero  for  t  >,  0. 
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should  be  zero  for  all  t  »  0.  Thus  we  must  have 


X(t)  =  e 


=  0 


-at 


t  >  0 


t  =  0 


( B— 16) 


That  is  we  must  understand  the  series  C B— 14)  to  mean  the  function  (B-16). 

"  (s) 

Similarly  any  function  }  b  6V  *  (t)  must  be  understood  as  identical  with  any 

s=0  s  r  s 

function  whose  Laplace  transform  is  \  b  p  (for  p  *  0)*  except  possibly  at 

s=0  s 

t  =  0.  Now  we  find  that,  using  (B-8) 


.nP(s)^x  _  ,  st\  s 

C  6  (t)  “  (  )  Ts~~V. 


=  0 


(0  <  n  s) 


.  (B-I7) 


n  > 


and  so 


OO  00 

[■  y » i(,)(t) .  <-)n  y  b 

/  i  s  /  i  n !  s+n 

s=0  s=0 


(t) 


(n  -  0  +  -)  .  (B-l  8) 


Consequently 


i.’t 


t-X» 


l  bs5<3><t) 
s»0 


0  for  any  positive  or  zero  n 


(B-l 9) 


We  conclude  therefore  that  any  infinite  series  of  6  functions  is  a  representa¬ 
tion  of  a  function  which  is  transcendentally  small  (t  -*■  °°)  compared  with  the 


Note  the  requirement  p  >  C 
all  t  >  0  ,  including  t  ■ 
p  ■  0  as  well  as  at  other  values  of 
at  any  p  including  p  ■  0 


If  we  want  the  function  to  be  equivalent  for 
,  then  their  transforms  must  be  identical  at 


p  .  The  transform  of  6^s^(t) 
and  from  the  theorem  of  equation  (5) 


is 


£St 

p+0 


y(p) 


OO 

/ 


y(t)dt 


If  the  M(a)  was  zero  in  the  above  example  the  transform  of 
p  ■  0  would  not  be  the  same  as  that  of  X(t)  . 


-at 


at 
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negative  powers  of  t  .  Thus  it  can,  for  example,  be  a  function  which  is  zero 
for  t  greater  than  some  finite  value,  or  a  function  which  behaves  like 
exp|-at^|  (a  >  0,  b  >  0)  as  t  tends  to  infinity.  The  latter  possibility  is 
illustrated  by  the  three  representative  cases  b  *  I,  2  for  which  we  have  the 
following  transforms: 


i/’|exp(-at^)| 


'©{' '  er£(^r)} 


(B-20) 


S'W-at)  }  '  TFT75 


(B-21) 


*{exp(-at2)}  -  A  „„(£_)  j,  - 


(B-22) 


where  erf(x)  is  the  error  function*.  Of  these  we  cannot  identify  the  first 
with  a  power  series  in  p  for  if  we  continue  it  analytically  along  the  negative 
real  axis  it  becomes  complex.  The  same  is  true  for  any  value  of  b  less  than 
one  and  so  we  can  say  that  |exp  -at*3 }  (a  >  0,  1  >  b  >  0)  is  not  equivalent  to  a 
series  of  6  functions.  The  other  two  possibilities  are  quite  acceptable;  and 
so  finally  we  can  say  that  any  infinite  series  of  5  functions  is  a  representa¬ 
tion  of  a  function  which  either  decays,  as  t  tends  to  infinity,  as  a  simple 

-gtl  ,  ,  , 

exponential  e  (or  more  rapidly),  or  is  zero  for  t  greater  than  some  finite 
value.  The  function  X(t)  of  equations  (B-J4)  and  (B-I6)  is  an  example  of  the 
first  possibility.  An  example  of  the  second  is  the  following. 

It  is  well  known  that  the  inverse  transform  of  e  pIq(p)  is 


*4(2  -  t) 


H(2  -  t). 


(B-23) 


However,  using  the  power  series  expansions  of  e  p  and  Ig(p)  we  find  that  it 


2  / 

Note  that  eX  jl  -  erf(x)}  ~  <1 - -j  +  ...  I  and 

4x  /  2x  ) 


erf  (x)  -  -^=  |  x  -  +  . . . 
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«(t)  -  <S0)(t)  +  |  6 (2)  (t)  -  ~2  6(3)(t)  +  6(4)(t)  -  ...  .  (B-24) 

The  difference  between  these  two  functions  must  be  a  null  function  and  so  inte¬ 
grating  both  as  before  between  0  and  t  we  find  that 

H(t)  -  6  (t)  +  ±  60)(t)  -  £  6(2)(t)  +  rA  6 (3)  (t)  -  ... 

=  g (t)H(2  -  t)  +  g(2) | 1  -  H(2  -  t)|  (B-25) 

where*  g(t)  =  ^sin  *(I  -  t)  -  .  (B-26) 

Thus 

C(t)  =  -  6(t)  +  |  6 (0  (t)  -  £  6 (2)  (t)  +  <S(3)(t)  -  ... 

=  g(t)H(2  -  t)  +  1  -  H(2  -  t)  -  H(t)  .  (B-27) 

This  infinite  series  of  <5  functions  ?(t)  therefore  represents  the  function  on 
the  right  hand  side  of  (B-27)  which  is  zero  for  t  *  0  and  for  t  >  2  .  More 

generally,  since  6  (•£■)  =  k6(t)  ,  we  have 


-  kfi(t)  +  |  k2<$0)(t)  -  f2kV2)(t)  +  ^kV3)(t)  -  ... 

g(£)H(2k  -  t)  +  1  -  H(2k  -  t)  -  H(t)  . 


(B— 28) 


It  may  well  be  that,  for  some  given  series  of  6  functions,  a  combination 
of  a  finite  number  of  series  such  as  (B— 14)  and  (B-28)  can  be  found  which  only 
differs  from  the  given  series  in  a  finite  number  of  terms.  In  such  a  case  there 
is  no  point  in  finding  the  representation  of  this  remainder  finite  series  as  a 
linear  combination  of  an  infinite  basic  set  of  infinite  series. 


The  value  of  sin  '(1  -  t) 


is  that  in  the  range  j  4  sin  *(1  -  t)  s 


J 
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N>|  X 


Real  negative 

P 

Real  positive 

P 

Purely  imaginary 

P 

X 

2 

X 

1 

M 

■n 

Z(±ix) 

0 

1 

1 

1 

0.1 

1.1241  *  0.4753i 

0.8024 

0.8319  *  0. 1723i 

0.2 

0.7677  T  0.7371i 

0.7315 

0.7276  ¥  0. 1886i 

0.3 

0.4605  ¥  0.6368i 

0.6901 

0.6650  *  0. 1 793i 

0.4 

0.3240  ¥  0.4781i 

0.6622 

0.6250  +  0. 1 650i 

0.5 

0.2768  ¥  0.3536i 

0.6418 

0.5979  *  0. 1507i 

0.6 

0.2672  *  0. 2647i 

0.6263 

0.5788  ¥  0. 1 378i 

0.7 

0.2731  ¥  0.2012i 

0.6139 

0.5648  ¥  0. 1 264i 

0.8 

0.2798  ¥  0. 1 550i 

0.6039 

0.5541  T  0. 1 1 65i 

0.9 

0.2959  ¥  0. 1 207i 

0.5955 

0.5459  *  0. 1078i 

1.0 

0.3141  ¥  0.09482i 

0.5884 

0.5394  +  0. 1003i 

2 

0.4110  *  0.01063i 

0.5512 

0.5130  ¥  0.05769i 

3 

0.4475  ¥  0. 001360i 

0.5362 

0.5063  ¥  0.04000i 

4 

0.4634  ¥  0.0001 794i 

0.5280 

0.5037  ¥  0.03050i 

5 

0.4718  ¥  0.00002393i 

0.5229 

0.5024  ¥  0.02460i 

10 

0.4868  ¥  0.0 . i 

0.5199 

0.5006  ¥  0.01 245i 

0.5 


Table  2 

VALUES  OF  pr  AND  T (r) ( 1 ) /r (r  +  ! ) 


(see  Appendix  A) 


■ 

Pr 

r(r)(i)/r(r  +  i) 

1 

-0.577215... 

-0.577215... 

2 

1.644934... 

0.989056... 

3 

-1.202057... 

-0.907479... 

4 

1.082323.. . 

0.981728... 

5 

-1.036928... 

-0.981995... 

6 

1.017343... 

0.993149... 

7 

-1.008349... 

-0.996002... 

8 

1.004077... 

0.998106... 

9 

-1.002008. . . 

-0.999025... 

10 

1.000995... 

0.999516... 

1 1 

-1.000494... 

-0.999757... 

12 

1.000246... 

0.999878... 

13 

-1.000123... 

-0.999939... 

14 

1.000061... 

0.999970... 

15 

-1.000031... 

-0.999985... 

16 

1.000015... 

0.999992.. . 

17 

-1.000008... 

-0.999996... 

18 

1.000004... 

0.999998... 

19 

-1 .000002. .. 

-0.999999... 

20 

1.000001... 

1.000000... 

<->*  £  ±  (r  » 

S=1 


p 


Y 


where  y  is  Euler's  constant 


(see  Appendix  A) 


LIST  OF  SYMBOLS 


A 


D 


inertia  matrix  (structural  +  aerodynamic) 
aerodynamic  inertia  matrix 
see  equation  (2) 
aerodynamic  damping  matrix 
=  B(-) 

aerodynamic  stiffness  matrix 
=  C(0) 

=  C(-) 

structural  damping  matrix 

_  _d_ 
dt 


E 

G(vt) 

G(p) 

G(vt) 


H(t) 


I 


Vz) 

K(t) 


Kq(2),  Kj (z) 

Ko 

K 

r 


S 

S(X) 


structural  stiffness  matrix 

constituent  of  indicial  aerodynamic  matrix  -  see  equation  (2) 
Laplace  transform  of  G(t) 

the  difference  between  G(vt)  and  its  asymptotic  (t  ■+■  °°)  power 
scries  expansion 

right  hand  Heaviside  step  function  -  see  equation  (B-3) 
unit  matrix 

modified  Bessel  function  of  the  first  kind  of  order  zero 
indicial  aerodynamic  matrix  -  see  equations  (1)  and  (2) 
modified  Bessel  functions  of  the  second  kind 
see  equation  (2) ,  ■  - 

coefficient  matrices  in  the  Richardson*  approximation  to  G(vt) 
coefficient  matrices  in  the  expansion  of  G(p)  -  see  equation  (11) 
the  characteristic  matrix  -  see  equation  (33) 

coefficient  matrices  in  the  expansion  of  G(p)  -  see  equation  (11) 
see  equation  (12) 
see  equation  (12) 

■  X  *  -  I  (Appendix  A) 
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Sr,(v> 

S*aM 


T  (v) 
rs 

T*  (v) 
rs 

W. 

1 

w* 

1 

X(X) 


x0(t) 

VT) 


q(t) 

q 


8rs<X> 


t 

v 


wr(X) 


x£j(X) 


LIST  OF  SYMBOLS  (continued) 

coefficient  matrices  in  S/  *  (£)|  ~  see  equation  (39) 
see  equations  (50)  and  (51) 
see  equation  (12) 

coefficient  matrices  in  M  1 (£•)  -  see  equation  (35) 
see  equation  (50) 

=  ( adjoint  of  M(u) *  j  -  (see  equation  (34)) 

'  Ai=X./u 

see  equations  (48)  and  (49)  1 

the  infinite  lower  triangular  matrix  with  non-zero  elements  given 
by  equation  (A-22) 


=  £S 


dJx, 


the  fundamental  solution 


dT 


J 


real  constant  column  matrices  which  are  the  coefficients  of 
6^(t)  in  an  instantaneous  excitation 

Laplace  transform  variable 

constant  in  Richardson*  approximation  to  G(vt) 
column  matrix  of  generalised  coordinate 

semi-amplitude  of  q  if  motion  is  simple  harmonic  of  infinite 
duration 

elements  of  the  lower  triangular  matrix  S  -  see  Appendix  A 
time 

airspeed  divided  by  a  reference  speed 

elements  in  the  first  column  of  I  +  S,  see  equations  (A-26)  to 
(A-28) 

elements  of  X(X) 


r(z) 

S(p) 


the  Ganna  Function 

■  K,(f]/{K0(f)  +  the  Theodorsen  function 
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a! 

l 


(rs) 


Y 

6(0 
6  (n) 

c(0 


(t) 


T 

X 


x<t) 


^(z) 


*''l 

f(p) 


^'W>| 


LIST  OF  SYMBOLS  (continued) 

left  hand  characteristic  vector 

see  equation  (43) 

right  hand  characteristic  vector 

*  0,577215...  ,  Euler's  constant 

right  hand  Dirac  delta  function,  see  Appendix  B 

dn6 

dtn 

see  equation  (B-23) 
see  equation  (29) 
see  equation  (29) 


0  (t) 
m 

=  &  l(logmp) 

X. 

1 

characteristic  value 

Un 

see  equations  (A-44)  and  (A- 

V 

=  w/v,  frequency  parameter 

cn(x) 

-  jj(D  +  X)njr(z)|j 

(p. 

=  -  Y 

no 

|pr 

-  <->r  I  ~  (r  >,  2) 
s=  1  sr 

al 

*  X  -  Y 

a 

r 

=  Pr  (r  *■  2) 

time  multiplied  by  (reference  speed/reference  length) 

arbm^rary  constant  used  in  getting  inverse  Laplace  transform  of 
logF*  see  Appendix  A 

see  equation  ( B— 14) 
dr  ,  , 

*  ~dz*  108  ’  the  functions 

frequency  multiplied  by  (reference  length/reference  speed) 

*  Laplace  transform  of  f(t) 

■  inverse  Laplace  transform  of  f(p) 

-a'jf(t)} 
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LIST  OF  SYMBOLS  (concluded) 

=  n !  /  r  •'  (n  -  r)  ? 

dnf 

dx 


signifies  the  imaginary  part  of  x 
signifies  the  real  part  of  x 
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